Single-cell 'omic profiles of human aortic endothelial cells in vitro and human atherosclerotic lesions ex vivo reveal heterogeneity of endothelial subtype and response to activating perturbations.

Heterogeneity in endothelial cell (EC) sub-phenotypes is becoming increasingly appreciated in atherosclerosis progression. Still, studies quantifying EC heterogeneity across whole transcriptomes and epigenomes in both in vitro and in vivo models are lacking. Multiomic profiling concurrently measuring transcriptomes and accessible chromatin in the same single cells was performed on six distinct primary cultures of human aortic ECs (HAECs) exposed to activating environments characteristic of the atherosclerotic microenvironment in vitro. Meta-analysis of single-cell transcriptomes across 17 human ex vivo arterial specimens was performed and two computational approaches quantitatively evaluated the similarity in molecular profiles between heterogeneous in vitro and ex vivo cell profiles. HAEC cultures were reproducibly populated by four major clusters with distinct pathway enrichment profiles and modest heterogeneous responses: EC1-angiogenic, EC2-proliferative, EC3-activated/mesenchymal-like, and EC4-mesenchymal. Quantitative comparisons between in vitro and ex vivo transcriptomes confirmed EC1 and EC2 as most canonically EC-like, and EC4 as most mesenchymal with minimal effects elicited by siERG and IL1B. Lastly, accessible chromatin regions unique to EC2 and EC4 were most enriched for coronary artery disease (CAD)-associated single-nucleotide polymorphisms from Genome Wide Association Studies (GWAS), suggesting that these cell phenotypes harbor CAD-modulating mechanisms. Primary EC cultures contain markedly heterogeneous cell subtypes defined by their molecular profiles. Surprisingly, the perturbations used here only modestly shifted cells between subpopulations, suggesting relatively stable molecular phenotypes in culture. Identifying consistently heterogeneous EC subpopulations between in vitro and ex vivo models should pave the way for improving in vitro systems while enabling the mechanisms governing heterogeneous cell state decisions.

Single cell 'omic profiles of human aortic endothelial cells in vitro and human atherosclerotic lesions ex vivo reveals heterogeneity of endothelial subtype and response to activating perturbations.

Objective: Endothelial cells (ECs), macrophages, and vascular smooth muscle cells (VSMCs) are major cell types in atherosclerosis progression, and heterogeneity in EC sub-phenotypes are becoming increasingly appreciated. Still, studies quantifying EC heterogeneity across whole transcriptomes and epigenomes in both in vitro and in vivo models are lacking. Approach and Results: To create an in vitro dataset to study human EC heterogeneity, multiomic profiling concurrently measuring transcriptomes and accessible chromatin in the same single cells was performed on six distinct primary cultures of human aortic ECs (HAECs). To model pro-inflammatory and activating environments characteristic of the atherosclerotic microenvironment in vitro, HAECs from at least three donors were exposed to three distinct perturbations with their respective controls: transforming growth factor beta-2 (TGFB2), interleukin-1 beta (IL1B), and siRNA-mediated knock-down of the endothelial transcription factor ERG (siERG). To form a comprehensive in vivo/ex vivo dataset of human atherosclerotic cell types, meta-analysis of single cell transcriptomes across 17 human arterial specimens was performed. Two computational approaches quantitatively evaluated the similarity in molecular profiles between heterogeneous in vitro and in vivo cell profiles. HAEC cultures were reproducibly populated by 4 major clusters with distinct pathway enrichment profiles: EC1-angiogenic, EC2-proliferative, EC3-activated/mesenchymal-like, and EC4-mesenchymal. Exposure to siERG, IL1B or TGFB2 elicited mostly distinct transcriptional and accessible chromatin responses. EC1 and EC2, the most canonically 'healthy' EC populations, were affected predominantly by siERG; the activated cluster EC3 was most responsive to IL1B; and the mesenchymal population EC4 was most affected by TGFB2. Quantitative comparisons between in vitro and in vivo transcriptomes confirmed EC1 and EC2 as most canonically EC-like, and EC4 as most mesenchymal with minimal effects elicited by siERG and IL1B. Lastly, accessible chromatin regions unique to EC2 and EC4 were most enriched for coronary artery disease (CAD)-associated SNPs from GWAS, suggesting these cell phenotypes harbor CAD-modulating mechanisms. Conclusion: Primary EC cultures contain markedly heterogeneous cell subtypes defined by their molecular profiles. Surprisingly, the perturbations used here, which have been reported by others to be involved in the pathogenesis of atherosclerosis as well as induce endothelial-to-mesenchymal transition (EndMT), only modestly shifted cells between subpopulations, suggesting relatively stable molecular phenotypes in culture. Identifying consistently heterogeneous EC subpopulations between in vitro and in vivo models should pave the way for improving in vitro systems while enabling the mechanisms governing heterogeneous cell state decisions.

Meta-Analysis of Smooth Muscle Lineage Transcriptomes in Atherosclerosis and Their Relationships to In Vitro Models.

BACKGROUND: Vascular smooth muscle cells (VSMC) exhibit phenotypic plasticity in atherosclerotic plaques, and among other approaches, has been modeled in vitro by cholesterol loading. METHODS: Meta-analysis of scRNA-seq data from VSMC lineage traced cells across five experiments of murine atherosclerosis was performed. In vivo expression profiles were compared to three in vitro datasets of VSMCs loaded with cholesterol and three datasets of polarized macrophages. RESULTS: We identified 24 cell clusters in the meta-analysis of single cells from mouse atherosclerotic lesions with notable heterogeneity across studies, especially for macrophage populations. Trajectory analysis of VSMC lineage positive cells revealed several possible paths of state transitions with one traversing from contractile VSMC to macrophages by way of a proliferative cell cluster. Transcriptome comparisons between in vivo and in vitro states underscored that data from three in vitro cholesterol-treated VSMC experiments did not mirror cell state transitions observed in vivo. However, all in vitro macrophage profiles analyzed (M1, M2, and oxLDL) were more similar to in vivo profiles of macrophages than in vitro VSMCs were to in vivo profiles of VSMCs. oxLDL loaded macrophages showed the most similarity to in vivo states. In contrast to the in vitro data, comparison between mouse and human in vivo data showed many similarities. CONCLUSIONS: Identification of the sources of variation across single cell datasets in atherosclerosis will be an important step towards understanding VSMC fate transitions in vivo. Also, we conclude that cholesterol-loading in vitro is insufficient to model the VSMC cell state transitions observed in vivo, which underscores the need to develop better cell models. Mouse models, however, appear to reproduce a number of the features of VSMCs in human plaques.

Transcriptional Profiling Suggests T Cells Cluster around Neurons Injected with Toxoplasma gondii Proteins.

Toxoplasma gondii's tropism for and persistence in the central nervous system (CNS) underlies the symptomatic disease that T. gondii causes in humans. Our recent work has shown that neurons are the primary CNS cell with which Toxoplasma interacts and which it infects in vivo This predilection for neurons suggests that T. gondii's persistence in the CNS depends specifically upon parasite manipulation of the host neurons. Yet, most work on T. gondii-host cell interactions has been done in vitro and in nonneuronal cells. We address this gap by utilizing our T. gondii-Cre system that allows permanent marking and tracking of neurons injected with parasite effector proteins in vivo Using laser capture microdissection (LCM) and RNA sequencing using RNA-seq, we isolated and transcriptionally profiled T. gondii-injected neurons (TINs), Bystander neurons (nearby non-T. gondii-injected neurons), and neurons from uninfected mice (controls). These profiles show that TIN transcriptomes significantly differ from the transcriptomes of Bystander and control neurons and that much of this difference is driven by increased levels of transcripts from immune cells, especially CD8+ T cells and monocytes. These data suggest that when we used LCM to isolate neurons from infected mice, we also picked up fragments of CD8+ T cells and monocytes clustering in extreme proximity around TINs and, to a lesser extent, Bystander neurons. In addition, we found that T. gondii transcripts were primarily found in the TIN transcriptome, not in the Bystander transcriptome. Collectively, these data suggest that, contrary to common perception, neurons that directly interact with or harbor parasites can be recognized by CD8+ T cells.IMPORTANCE Like other persistent intracellular pathogens, Toxoplasma gondii, a protozoan parasite, has evolved to evade the immune system and establish a chronic infection in specific cells and organs, including neurons in the CNS. Understanding T. gondii's persistence in neurons holds the potential to identify novel, curative drug targets. The work presented here offers new insights into the neuron-T. gondii interaction in vivo By transcriptionally profiling neurons manipulated by T. gondii, we unexpectedly revealed that immune cells, and specifically CD8+ T cells, appear to cluster around these neurons, suggesting that CD8+ T cells specifically recognize parasite-manipulated neurons. Such a possibility supports evidence from other labs that questions the long-standing dogma that neurons are often persistently infected because they are not directly recognized by immune cells such as CD8+ T cells. Collectively, these data suggest we reconsider the broader role of neurons in the context of infection and neuroinflammation.

Systems Genetics in Human Endothelial Cells Identifies Non-coding Variants Modifying Enhancers, Expression, and Complex Disease Traits.

The identification of causal variants and mechanisms underlying complex disease traits in humans is important for the progress of human disease genetics; this requires finding strategies to detect functional regulatory variants in disease-relevant cell types. To achieve this, we collected genetic and transcriptomic data from the aortic endothelial cells of up to 157 donors and four epigenomic phenotypes in up to 44 human donors representing individuals of both sexes and three major ancestries. We found thousands of expression quantitative trait loci (eQTLs) at all ranges of effect sizes not detected by the Gene-Tissue Expression Project (GTEx) in human tissues, showing that novel biological relationships unique to endothelial cells (ECs) are enriched in this dataset. Epigenetic profiling enabled discovery of over 3,000 regulatory elements whose activity is modulated by genetic variants that most frequently mutated ETS, AP-1, and NF-kB binding motifs, implicating these motifs as governors of EC regulation. Using CRISPR interference (CRISPRi), allele-specific reporter assays, and chromatin conformation capture, we validated candidate enhancer variants located up to 750 kb from their target genes, VEGFC, FGD6, and KIF26B. Regulatory SNPs identified were enriched in coronary artery disease (CAD) loci, and this result has specific implications for PECAM-1, FES, and AXL. We also found significant roles for EC regulatory variants in modifying the traits pulse pressure, blood protein levels, and monocyte count. Lastly, we present two unlinked SNPs in the promoter of MFAP2 that exhibit pleiotropic effects on human disease traits. Together, this supports the possibility that genetic predisposition for complex disease is manifested through the endothelium.

Sibling Rivalry in Myxococcus xanthus Is Mediated by Kin Recognition and a Polyploid Prophage.

UNLABELLED: Myxobacteria form complex social communities that elicit multicellular behaviors. One such behavior is kin recognition, in which cells identify siblings via their polymorphic TraA cell surface receptor, to transiently fuse outer membranes and exchange their contents. In addition, outer membrane exchange (OME) regulates behaviors, such as inhibition of wild-type Myxococcus xanthus (DK1622) from swarming. Here we monitored the fate of motile cells and surprisingly found they were killed by nonmotile siblings. The kill phenotype required OME (i.e., was TraA dependent). The genetic basis of killing was traced to ancestral strains used to construct DK1622. Specifically, the kill phenotype mapped to a large "polyploid prophage," Mx alpha. Sensitive strains contained a 200-kb deletion that removed two of three Mx alpha units. To explain these results, we suggest that Mx alpha expresses a toxin-antitoxin cassette that uses the OME machinery of M. xanthus to transfer a toxin that makes the population "addicted" to Mx alpha. Thus, siblings that lost Mx alpha units (no immunity) are killed by cells that harbor the element. To test this, an Mx alpha-harboring laboratory strain was engineered (by traA allele swap) to recognize a closely related species, Myxococcus fulvus. As a result, M. fulvus, which lacks Mx alpha, was killed. These TraA-mediated antagonisms provide an explanation for how kin recognition specificity might have evolved in myxobacteria. That is, recognition specificity is determined by polymorphisms in traA, which we hypothesize were selected for because OME with non-kin leads to lethal outcomes. IMPORTANCE: The transition from single cell to multicellular life is considered a major evolutionary event. Myxobacteria have successfully made this transition. For example, in response to starvation, individual cells aggregate into multicellular fruiting bodies wherein cells differentiate into spores. To build fruits, cells need to recognize their siblings, and in part, this is mediated by the TraA cell surface receptor. Surprisingly, we report that TraA recognition can also involve sibling killing. We show that killing originates from a prophage-like element that has apparently hijacked the TraA system to deliver a toxin to kin. We hypothesize that this killing system has imposed selective pressures on kin recognition, which in turn has resulted in TraA polymorphisms and hence many different recognition groups.